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Abstract 

This paper investigates existence of the nonstandard Pade approximants introduced by 
Cherkaev and Zhang in |2| for approximating the spectral function of composites from ef- 
fective properties at different frequencies. The spectral functions contain microstructure infor- 
mation. Since this reconstruction problem is ill-posed [9], the well-performed Pade approach 
is noteworthy and deserves further investigations. In this paper, we validate the assumption that 
the effective dielectric component of interest of all binary composites can be approximated by 
Pade approximants whose denominator has nonzero power one term. We refer to this as the 
nonstandard Pade approximant, in contrast to the standard approximants whose denominator 
have nonzero constant terms. For composites whose spectral function assumes infinitely many 
different values such as the checkerboard microstructure, the proof is carried by using clas- 
sic results for Stieltjes functions. For those with spectral functions having only finitely many 
different values, we prove the results by utilizing a special product decomposition of the coef- 
ficient matrix of the Pade system. The results in this paper can be considered as an extension 
of the Pade theory for Stieltjes functions whose spectral function take infinitely many different 
values to those taking only finitely many values. In the literature, the latter is usually excluded 
from the definition of Stieltjes functions because they correspond to rational functions, hence 
convergence of their Pade approximants is trivial. However, from an inverse problem point 
of view, our main concern is the existence of the nonstandard Pade approximants, rather than 
their convergence. The results in this paper provide a mathematical foundation for applying the 
Pade approach for reconstructing the spectral functions of composites whose microstructure is 
not a priori known. 
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dard Pade approximants; Stieltjes functions; Inverse homogenization. 
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1 Introduction 



In the seminal paper by Golden and Papanicolaou[l], components of the effective dielectric matrix 
of a two-pahse composite with statistically homogeneous microstructure and isotropic constituents 
are represented as Stieltjes integrals with positive Borel measure supported in [0, 1]. We refer to 
this integral representation formula as IRF throughout the paper. In the IRF, the contrast of the 
constituents, represented as the ratio (denoted by h) of the two dielectric constants, stays in the 
integrand while the information of microstructure is contained in the positive Borel measure. To 
be more specific, it is shown in [ 1 ] that for s := outside [0, 1] on the complex plane, the (i, k)- 
component of the effective dielectric tensor, i, k = 1,2, 3 of a dielectric composite, denoted by e* k 
can be represented as 

e ik = e l (tik ~ I —^dfiik(z) 
V Jo s-z 

where h := |^-, is Kronecker delta, and d^ik(z) is a positive Borel measure. Throughout the 
paper, we adapt the notation := and =: such that 

A := B o A is defined by B 
A=: B o Bis defined by A 

Following the convention in literature, we introduce F(s) function 

e* f 1 1 

Fik(s) := <5j fc - — = / dfi ik (z) (1) 

ei Jo s- z 

Note that h is complex-valued for lossy materials. Applying geometric series expansion of 
around s = oo (corresponding to a homogeneous material), ([I]) becomes 

X 

Fik(s) = ^ 
m=0 

with c m being the m-th moment of the measure, i.e. c m := z m fiik(dz). 

It is also shown in 0] that c m , m = 0, 1, 2, • • • , is related to the (m + 1) -point correlation 
function of the characteristic function of the region occupied by one of the constituents. This link 
has been exploited for estimating the bounds on effective dielectric tensor by using information on 
m-point correlation functions such as volume fraction (m = 1) and two-point correlation function. 
This approach of estimating bounds relies on the observation from ([!]) that the i^fc(s) is an analytic 
function for s € C \ [0, 1], where C is the complex plane. This provides a mathematical scheme 
for retrieving the moments c m when the dielectric constant of the one of the constituents varies 
with frequency. In 0, it was demonstrated analytically and numerically that the construction of 
moments of the measure from data of e* k at various frequencies is a stable process, unlike the ill- 
posed process of reconstructing the measure itself. In that paper, the reconstruction of moments 
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was carried out by solving a truncated linear system if all the data points are outside the unit 
circle and by a Gaussian quadrature method. In all the tested numerical examples, the first 5-10 
moments can be reconstructed with high accuracy by both methods. These moments contain the 
microstructural information and the resulting truncated series also provide a very accurate way for 
calculating the effective dielectric constant for all frequencies(extrapolation) [3]. In [2], instead of 
reconstructing the coefficients of the truncated Stieltjes series, the coefficients of the nonstandard 
Pade approximants of the series in ([I]) was reconstructed from data of e* k at different frequencies. 
This method is able to reconstruct the first 10-15 moments with high accuracy and the resulting 
Pade appriximant provides a very accurate way for calculating the effective dielectric constant for 
all frequencies(extrapolation). It is worth noticing that all the real- valued poles of the resulting Pade 
approximant lie along the graph of the true spectral functions in all the tested numerical examples. 
To be more specific, in 0, Pade approximants for the Stieltjes integral in ([T]) with the following 
form: 

was used for recovering the moments. They assume F(s) has at least one pole and normalize b\ 
to be 1, unlike the standard Pade approximants for which bo is normalized to 1. We refer to (|2]) as 
nonstandard Pade apprximants. The main reason prompting the choice of b\ ^ is beecause 
for a checkerboard microstructure, the spectral function has a pole at s = 0. For ease of notation, 
we will replace e* k with e* hereafter. \i is assumed to take at least two different values because 
otherwise F will be identically zero as implied by ([TJ. For N different frequencies oj\, ■ ■ ■ ,ojn, 
there are N data pairs (sk, dk), where Sk = i- €2 ^y £1 ( Uk ) an d dk = e*(sfc), k = 1, ■ • ■ , N. The 

• , dp, bo, »2 

ap + a lSk + a 2 s 2 k + ■■■ + a p s p k 

d k-—r—, , , 2 i ., q ,p + g + l<iV (3) 

bo + s k + b 2 s% H h b q s\ 

was solved by formulating it as a least-square minimization problem subject to the constraints 

0<A n <l, 0<a„<l, 0<^A 1 <1 (4) 

n 

for the residues A n and poles s n defined by the partial fraction decomposition of the nonstandard 
Pade approximant 

ao + a\s + 02S 2 + • • • + a p s p _ ^ A n 



following linear system of equations for ao, • • • , a p , 6q, 62, • • • ,b, 



b + 61 s + b 2 s 2 + • • • + b Q si , s - s n 

The moments were then computed by applying partial fraction decomposition to the Pade approx- 
imants and by correctly combining the poles and residues. See (21 |4j for details. 

In (21 @ Q , the fact that a Stieltjes function can be bounded up and below by standard Pade 
approximants was used for deriving bounds on the effective dielectric constants. 
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The main purpose of this paper is twofold. The first is to show that F(s) can always be approx- 
imated by a Pade approximant of the form in ^ for s G C \ [0, 1], i.e. to show that we can always 
normalize 61 to be 1 and to show this special form of Pade approximant has the accuracy-through- 
order property that will be defined in the next section. The second is to justify the constraints shown 
in @. 

This paper is organized as follows. We first define an auxiliary function /(£) and prove that 
it has accuracy-through-order standard Pade approximants in Section |2.1| In the same section, 
we also prove properties of the poles of the Pade approximants. Using these results, we show in 
Section [5] that F(s) always has accurate-through-order Pade approximant of the form as in §2§ for 
appropriate p and q. The major technique in the proof is fractorizing the coefficient matrix of the 
system of Pade equations into products of Vandermonde matrices and a diagonal matrix. In Section 
[4| we discuss the relevance of this work in materials science for binary composites. 

2 Existence of Pade approximants 

Since the representation formula for F(s) is a Stieltjes integral with non-decreasing distribution 
function fi and is analytic for s € C \ [0, 1], we may identify it with a Stieltjes function /(£) (8l by 
rewriting F(s) in terms of a new variable £ := — 1/s 

F(s) =: G(0 = f 1 ^ = : • /(0 < 5 > 
Jo 1 + Z K 

We refer to /(£) as the auxiliary function and would like to remark that the function can take 
finitely many or infinitely many different values, i.e. the definition of Stieltjes function we use here 
does not include the restriction on the number of different values n(z) takes. Note that the Stieltjes 
function /(£) is analytic for £ € C \ (— 00, —1] and has a power series expansion valid in |£| < 1, 
which reads 

00 

/(0=E^n("O ft - (6) 
n=0 

For later use, we also include the expansion of G in the same region on the £-plane 

00 

G(0 = J2^(~0 n+1 - (7) 

n=0 

Here we use — £ because it is easier to work with the positive coefficients which can be easily 
verified as the moments of the Borel measure /iin(B]),n = 0,l,2,---. 
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2.1 Existence of standard Pade approximants for the auxiliary function /(£) 



Recall the standard [L/M] Pade approximant for /(£) in ^ is a rational function of the following 
form 

[L/M]{ -° - b + bl (-o + ... + bM (-OM 

such that 60 can be normalized to be 1 and the rest of the coefficients satisfy the system of Pade 
equations (H 



(8) 



HL-M+1 


HL-M+2 ■ 






bu 






HL-M+2 


HL-M+3 ■ 


fJ-L+1 




&M-1 






/'L 




■ HL+M-1 




b l 




HL+M 



(9) 



ao = no 
oi = Mi + & i^o 

min(L,A/) 
k=l 

For consistency, we define /ij = for j < 0. The [L/M] Pade approximant is said to have the 
accuracy-through-order property if it satisfies 

f(0-[L/M) = O(e +M+1 ). (10) 

for |£| within radius of convergence. We would like to remark that for any given function f(z), 
polynomials 

p [L/M]^ q [L/M] 

(z) of degrees L, M, respectively, can always be found so that 



JL/M] 



(z)f{z)-p^ M \z) = 0{z 



L+M+l\ 



(11) 



However, it is well known that (111 does not necessarily imply p^ L / M \z) / q^ L ^ M \z) = f(z) + 
0{z L+M+l ) EQ. Following Baker, we say that the Pade approximant does not exist if the accuracy- 
through-order requirement is not satisfied. On the other hand, it can be shown that IH(p.7) the 
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p[L/M](t\ 

rational function murr, with p[ L / A/ ] and Q^ L / M ^ defined as 

Q[L/M](_£) 



ML-M+1 
ML-M+2 

AiL-1 

ML 



ML-M+2 
ML-M+3 

Ml 



L-AI 



L-M+l 



M+i-l 



i=0 



ML+1 
ML+2 

ML+M— 1 
L 

i=0 



g^(-e) 



ML-M+1 ML-M+2 

ML-M+2 ML-M+3 

Ml-i Ml 

ML ML+1 

(-6 M (-O^- 1 



ML 
ML+1 

ML+ Af-2 
ML+M-1 

(-0 



Ml+i 

ML+2 



ML+ Af-1 
ML+M 

1 



(13) 



is an accuracy-through-order Pade approximant of /(£) provided Q^^ M ^(0) 7^ 0. Note that 
Q[ L / M ](0)(Hankel determinant) is exactly the determinant of the coefficient matrix of ([8]>. Hence 
the condition Q[ L / M ](0) 7^ guarantees the accuracy-through-order property of the [L/M] ap- 
proximant [8](p.21). If q[ l / m 1(o) = 0, m can have infinitely many solutions or fail to have 
any solution, i.e. the [L/M] Pade approximant may not exist. We would like to remark that 
Q\ l / m \(Q^ = does not imply non-existence of accuracy-through-order Pade approximants. For 
example, for a rational function with degree p numerator and degree q denominator, it is known 
that Q[ L / M 1(0) = if L > p and M > q but obviously, the Pade approximants are the rational 
function itself for these cases. 

Because the numerical algorithm in (U is based on the reconstruction of the Pade approximant 
of F(s), rather than F(s) itself, it is necessary to prove the existence of Pade approximant for 
defined in (|5]). To achieve this, we first show the existence of standard Pade approximants for the 
Stieltjes function /(£) in the following theorem. 



Theorem 2.1. The [L/M] Pade approximant for /(£) exists for all L — M + 1 > 0, L ^ - 1 and 
M > such that b = 1. 

Proof. We first consider the case where the non-decreasing function p(z) in ^ takes infinitely 
many different values. Using Theorem 5.2.1 in [8], which states that Q^ L / M \0) 7^ for all L — 
M + 1 > and M > 1, the existence of [L/M] approximants follows directly. For M = 0, [L/0] 
clearly exists because [L/0] = ^^ = oMj( — ■ 
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For the case in which /j,(z) takes n + 1 different values for some non-negative finite integer 
n, we prove the theorem as follows. Let z\, 22, • • • z n be the n different points on [0, 1] such that 
< z\ < Z2 < 23 < • • • < z n and /jl(z) changes values only at these points, i.e. d/j,(z) is a finite 
sum of Dirac measure sitting at z% with strength Aj. 

Case 1: n = If n = 0, then /(£) = so [L/M] = for all L - M + 1 > and M > 0. 

Case 2: n = 1 If z\ = 0, i.e. /(£) takes two different values and the value change occurs at 
z\ by the amount of \\, then /(£) = Ai > 0. Clearly, the Pade approximant exits. If z\ > 0, 
/(£) = Tq A_ so [L/M] = Tq ^ ? for all L - M + 1 > and M > 1 by Theorem 1.4.4 in QQ and 

[L/0] clearly exists because [X/0] = X^=o ^j^ - J - 

Case 3: n > 2 Denote the amount of value change at Zi by Aj > 0, i = 1, 2, • • • , n. With this, 
/(£) can be written as 

i=i JS 

Suppose zi 7^ 0, i.e. Zj, i = 1 ~ n are all non-zero. Then /(£) is a rational function with degree 
n — 1 numerator and degree n denominator. By the characterization theorem of rational function 
(Theorem 1.4.4) in [8], we know that [L/M](-g) = /(£) if L > n - 1 and M > n, i.e. the 
Pade approximant of a rational function is itself if both L and M are large enough. Also, the 
constant term of the denominator is obviously non-zero. This takes care of those [L/M] such that 
L - M + 1 > and M > n. 

To complete the proof for z\ 7^ 0, we consider [L/M] such that L—M+l > and 1 < M < n. 
For this case, we exploit the the structure of moments /j,q = Y^=i \\ = Y17=i ^i z i' ^ = 
1, 2, 3 . . . to rewrite the coefficient matrix of ([8]) as 

r n 

i=i 

n 

E^ 

8=1 
n 

E XiZ i 
- i=l 
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-M+l \ „L-M+2 



1 E^ 

i=i 

n 

E^ 



-M+2 \ „£-M+3 



i=l 



E A ^ 

i=l 



8=1 



V T AV, 



(14) 



Em l+1 



E 

i=i 



A,;z, 



L+M-l 



where 



Z2 



y M-l 
J l 

JW-1 



V := 



1 Z M 



y M-l 



(15) 



A :-- 



\\Z 






L-M+l 
1 



\ 2 z 2 



L-M+l 







\ 3 Z 



M+l 



JVf-1 



^n z n 



nxM 



L-M+l 



(16) 



The superscript T denotes matrix transpose. Note that the null space Af(V) = E R because 
V is a rectangular Vandermonde matrix with re rows and M columns such that re > M. Let x 
be a vector in M M so that V T AVa; = 0. Consider the inner product (a;, y) := J2k=i x iVi f° r 
x = (xi, • • • ,x n ) T , y = (yi, ■ ■ ■ ,y n ) T , then 



= (V T AVx,x) 



(AVs, Va;) > min \z. 



L-M+l 



(V£c,Va;) 



Since min \;zf M+1 > 0, we have Yx 



0. Hence x = G R M . This shows that the 

coefficient matrix in ^ is non-singular so the [L/M] approximant exists for L — M + 1 > and 

1 < M < n if Zi 7^ for every i. 



If z\ = 0, /(^) is a rational function with degree n — 1 numerator and degree re— 1 denominator. 
The characterization theorem of rational function (Theorem 1.4.4) in [8] implies that [L/M] (— £) = 
/(£) if L > re — 1 and M > re — 1. This takes care of L — M + 1 > 1 and M > re — 1, represented 
by the triangular region marked by x in Table[T] For L — M+l = and M < re — 1, the coefficient 
matrix in ([8]) can be decomposed in the same fashion as in ( 14 1 with A replaced by 



A : = 



Ai 

A 2 

A 3 









A, 



(17) 



By an argument exactly the same as for the z\ ^ case, the existence of the [L/M] Pade approx- 
imant for L — M + 1 = and M < re — 1 is proved; they correspond to the line formed by o in 
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Table 1: Schematic description for case n > 2 and z\ = 



M\L 





1 


2 


77,-3 


n — 2 


71—1 


77 


n + 1 





• 


• 


• 


• 


• 


• 


• 


• 


1 





• 


• 


• 


• 


• 


• 


• 


2 




o 


• 


• 


• 


• 


• 


• 


3 






o 


• 


• 


• 


• 


• 


n-2 








o 


• 


• 


• 


• 


n — 1 










o 


X 


X 


X 


n 















X 


X 


n + 1 














o 


X 



Table [T] For L - M + 1 > and M < n - 1, the coefficient matrix of ^ is 



X-M+2 



8=2 
n 



t=2 



i=2 



i=2 



L-M+3 



i=2 



E^ L e^ 



i=2 



E* 

i=2 
i=2 



£ 

i=2 



L+M-l 



V T AV, 



(18) 



where 



V 



A := 



A 2 4" A/+1 



\ 3 zt M+1 



1 ^2 
1 z 3 Z- 



7 2 
2 



1 Z M Z 



M 



1 Zn Z n 














z 2 

„Af-l 



„M-1 



,Af-l 



(n-l)xM 



(n-l)x(n-l) 



(19) 



(20) 
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Since V is a Vandermonde matrix with n — 1 rows and M columns such that n — 1 > M and 
Zi 7^ Zj for i ^ j, the null space AA(V) = 06 R . Because all the entries in A are positive, the 
existence of the [L/M] Pade approximant can be concluded in the same way as for the previous 
case. The corresponding region in the Pade table is marked in Table [T]by •. □ 

We apply this result to show the existence of accuracy-through-order Pade approximants for 

G(0- 



2.2 Standard Pade approximants for G(£ ) 

Theorem 2.2. The [L/M] Pade approximant for exists for all I — M > and M > such 
that bo = 1. 

Proof For I > 1, let [I - 1/M]/ denote the [I - 1/M] Pade approximant of /(£) such that 
(I — 1) — M + 1 > 0, I > 1 and M > 0. The proof follows immediately from the definition of 



G, Theorem|2.1|and definition of accuracy-through-order property ( 10) 



G(0 " (~0[L ~ l/M] f = -e (/(0 — [L — 1/M]/) = 0(^+ M+1 ) for |^| < 1. 

That is, (-£)[! - 1/M]/ is the [I/M] Pade approximant of G(f) for I > 1 and I > M > 0. If 
1 = 0, then M = because of the condition I > M. The [0/0] approximant of G is because 
G(£) - = /U (-£) + /Ui(-C) 2 H = O(^ 0+0+1 ). The last equality is due to the positivity of 

Mo- □ 

The algorithm presented in [2] for reconstructing spectral function from known measured data 
e*(£) for £ not necessarily in |£| < 1. This can be justified by the two theorems listed below 
regarding properties of Pade approximants of Stieltjes series. 

Theorem 2.3. The poles of the Pade approximants of f, denoted by [I/M]/, I — M + 1 > and 
M > are simple poles lying on (— oo, — 1) with positive residues. 

Proof. For the case where the measure \i takes infinitely many different values, the proof can be 
found in [8 ], Theorem 5.2.1, p201. When n takes (re + 1) < oo different values, the proof given in 
1H needs to be modified as follows. The key step is to show the interlacing property of the zeros of 
A$(x) defined in O. 



If M = 0, the Pade approximant has no pole. For n = 0, 1, the proof is the same as that for 
Theorem 12. II 

For n > 2, we first consider the case when z\ ^ 0, i.e. no jump in values at z = hence 
/(£) is a rational function with numerator of degree n — 1 and denominator of degree n. Let 
J := I — M > —1, and M > 1. Following the proof of Theorem 5.2.1 in we define 
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A { Q J) (x) = 1; for M = 1, 2, 3 • • ■ , define 



A 



fj, 1+ j + XH2+J 
fJ, 2 +J + 



H2+J + XfJ-3+J 
+ XHi+j 



HM+J + XfJ,M+J+l 
HM+J+l + XflM+J+2 



HM+J + XfJ.M+J+1 HM+J+l + X/IM+J+2 ■ ■ ■ H2M+J-1 + X/I 2 M+J 



(21) 



Note that A^ } (x) = Q^ L / M \-x) by applying (3.2) on p. 16 of [|8]l to ([Ti]) and that A^ } (0) = 
for M > n and J > — 1. Hence the zeros of A^? correspond to the poles of Pade approximant 



[L/M]. Unlike the proof for n = oo in J8J, Theorem 5.2.1, we only need to consider J > — 1 and 
1 < M < n because for J > —1 and M > n, the Pade approximant for /(£) is itself the poles 
of /(£) with finite n are apparently all simple and located in (— oo, — 1). Next, we need to show 



that the functions a[ J \x), • • • , A,J 1 'j' 1 (x) have the property such that if A 



then A^^_ x (x) and A^ ; +1 (x) have opposite signs for all 1 < M < n — 1. This can be shown 
by using Sylvester's identity and Frobenius identity in exactly the same way as in the proof in [8], 
Theorem 5.2. 1. In order to show the interlacing property of the zeros of (x), we need to prove 
the following properties 



.(•0 



{J), 

M 



for some x, 



A$(0) > 
A^(-oo) = (-!)«■ oo 



forM = !,••• ,n 



(22) 



by first noting that A^ (0) is the determinant of the matrix V T AV in jl4| and 

(V T AVy, y) = (AVy, Vy) > min A^f" M+1 (Vy, Vy) > 0, for y G R M \ {0} 

So det(V T AV) > 0, hence A^(0) > for M = 0, 1, • • • , n. For A^(-oo), we decompose 
the matrix corresponding to A^ (x) as 



Hi +J + X/J.2+J 

fJ, 2+ J + XfJ, 3+ J 



fJ-2+J + Xfl 3+ J 
+ Xfi 4+ J 



V-M+J + X^LM+J+l MM+J+1 + X^lM+J+2 



HM+J + XflM+J+1 
fJ-M+J+1 + X^M+J+2 

^2M+J-1 + X/J.2M+J 



V T A Y V 



with 



A^ := 



X lZ ^- M+1 (l + z lX ) 

\ 2 z^- M+l (l + z 2 x) 








n 



Similar as before, for a fixed x, there exists A and B such that 

A(Vy, Vy) > (V T A x Vy,y) = (A x Vy, Vy) > B(Vy, Vy), for y e M M \ {0} 

where 

A:= max A,-^~ M+1 (1 + ^x), B := min A,-«f~ w+1 (l + Zjx) 
jr" = 1 , — ,n ' j=l,"' ,n J 

As x — > — oo, we have A — >• — oo and f? — > — oo. Hence A^(x)=det(V T A x V) — >• (— oo) M as 
x — )• — oo. This completes the proof of (22i for Z\ 7^ 0. 

For z\ = 0, we first consider that case when L — M + l = (i.e. J = — 1) and n — 1 > M > 1. 



We have A^ } (0) = det(V' J 'AV) with A defined in dl7k. So A^ ; (0) > for the same reason as 



in the Z\ 7^ case. As for A^ M '(x) for fixed x, we have '(x) = Ai + ^ A,(l + xzi) — >■ —00 
as x — > —00 and for M > 2 



i=2 



Ag } (x)=det 



Ai 



\ [ 








+ V T A X .V 



MxM 



with V being the one defined in ( 19 ) and 



A, 



X 2 (l + z 2 x) 

X 3 (l + z 3 x) 










Therefore, for any y = (yi, • ' ' > Um) T / 0, we have 

Ai • 



A n (l + z n x) 



(23) 



(n— l)x(n— 1) 



Aiy? + ^(Vy,Vy) > 





r 


V 


I. 









+ V T A X .V y, y > X iy f + B(Vy, Vy) 



MxM 



with ^4 := max Aj(l + Zjx) and B := min Aj(l + Zix); hence A^(x) — > (— oo) M as 



i=2,--- ,n 



00. For L — M + l > 1, A$(0) = det(V T AV) and A M '(x) = det(V T A x V) with V 



(J), 



defined in ( 19 1, A defined in §20y and A x defined in p3| ) and the same argument can be applied for 
proving (22 ). The rest of the proof is the same as that in [8 ], Theorem 5.2. 1. □ 

This theorem implies the following theorem: 
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Theorem 2.4 ([8], Theorem 5.4.4, p225). Any sequence of [L^/Mf.] Fade approximants, Lf, — 
Mk + 1 > 0, of a Stieltjes series /(£) convergent in |£| < 1 converges uniformly to /(£) in 
P + (A), where "D + (A) ij se£ o/a/Z points in |£| < i? which are at least A distance away from 
the cut (— oo, —1] /or aTry positive numbers R and A. 

It other words, even though the Pade approximants are constructed from the Stieltjes series 
convergent only in |£| < 1, their validity of approximating the Stieltjes function /(£), which is 
analytic everywhere outside (— oo, —1], extends far beyond |£| < 1. 

In 0, Cherkaev et. al. define Sn equivalence of dielectric composites with isotropic con- 
stituents. According to their definition, two microstructures are Sn -equivalent if the the first N 
moments (/io to of the spectral measures fi^ and fi^} in the IRF ([Ti are identical. In the 

following, we show that Pade approximants provide a characterization of Sn -equivalent structures 
in the following sense: 



Theorem 2.5. Let Gi(£) and G2(£) be two Stieltjes functions defined as in (|5]> by the effective 
complex dielectric constant of Structure- 1 and Structure-2, respectively. Denote the corresponding 
[L/M] Pade approximant by [L/M]g 1 and [L/M]a 2 - Then Structure-1 and Structure-2 are Sn 
equivalent, N > 1, if and only if [L/M] Gl = [L/M\g 2 for all L + M < N such that L > M > 0. 



Proof. Since the coefficients of (— n > 0,in the power series expansions of G\ and G2 ((VJ) 
are given by the n th moment of the corresponding measure, Sn -equivalence implies that the power 
series expansions near £ = of G\ and G2 have the same first N + 1 terms (the constant term in 
the series is for both G\ and G2 )■ Since the [L/M] Pade approximant is completely determined 
by coefficients up to {-£) L+M , clearly [L/M] Gl = [L/M] Ga if L + M < N and if they both 
exist. By Theorem [TT] we conclude that [L/M] Gl = [L/M] G2 for L + M < N and L > M > 0. 



Conversely, suppose [L/M] Gl (-£) = [L/M] G2 (-£) for L + M = N. Due to Theorem [23 
all the poles of Pade approximants for Stieltjes series convergent in |£| < 1 lie on the real axis in 
the interval — 00 < £ < — 1, the power series expansion of this rational function 



00 



[L/M] Gl (-0 = [L/M] G2 (-0 = Y J C^ M K-O n - 

n=0 

is valid for |£| < 1. By accuracy-through-order property of the Pade approximants, cif^ = p^-i 

and 4f^ M ' = for n = 1 ~ ^ with ixh 1 and ^ being the coefficients of G\ and G2, 

respectively, as defined in ([7]). Hence the theorem is proved. □ 

3 Existence of nonstandard Pade approximants for F(s) 

The main result of this paper is to show the existence of nonstandard Pade approximants in (|2]). 
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Theorem 3.1. For any Stieltjes function 

F{s)= I -^dfi(z) (24) 



with non-decreasing spectral function p,(z) taking either infinitely or finitely many different values 
for z G [0, 1], there exits nonstandard Pade approximant of the form 

F ( S )-TTT= — I , h , h 2 ^ TT M — >6l^0,M>l. (25) 

b(s) b + bis + b 2 s z + • • • + b M s M 

Proof. We first consider the case in which fj,(z) takes infinitely many different values on [0, 1]. Let 



L = M - 1 and M > 1. According to Theorem [271] and Theorem |23j the [M - 1/M] Pade 
approximant of /(£) exists with simple poles lying on £ G (— oo, —1) with positive residues. The 

interlacing-of-zero property of (x) guarantees the existence of M different poles of [M — 
1/M]. By Theorem 2.3 it can be further expressed in terms of the partial fraction decomposition 

[M ~ 1/M], = =J2j~' Pi e (-°°. » * Pi 

for i + j, Xi > 0. (26) 

This induces the [M/M] approximant of G(£): 

[M/M] G = -i ■ [L/M] f = 7 s — ^ Pi G (-oo, -1), Aj > 0. (27) 

i=i 4 Pi 

We rewrite the expression above in terms of variable s := ^ G [0, 1) 

Noting — G [0, 1) and F(s) = G(£), the partial fraction representation above implies there exists 



a rational function approximation of F(s) with the following form 

,. , a(s) _ ap + ais + a 2 s 2 H h cm-i^*" 1 



where 6i 7^ 0. This is because the coefficient of s in the denominator of (28 1 is J2i=i Tlj^i 
which is non-zero for M > 2. For M = 1, this is obviously true. 

Next, we consider the case for fi taking (n + 1) different values, 1 < n < 00 with jumps at 
1 > zi, . . . , z n > of magnitude Ai, . . . , A n > 0. Then 

n X 

F ( s ) = E — w 
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For M > n, the [M — 1/M] Pade approximant of F(s) is itself and apparently (29 1 has the form 
(25 ). For M < n, the proof is the same as in the case where \i takes infinitely many different values 
because of the validity of Theorem |2.1| and Theorem [23] for this case. □ 



4 Conclusion and discussion 



In this paper, we prove that a Stieltjes function of the form ( 24 ) can be approximated by [M — 1 / M] 
non-standard Pade approximant p5j ). Most of the work in this paper is dedicated to the case 
where the spectral measure \i takes finitely many different values in [0, 1] because the proof for 
H that takes infinitely many different values in [0, 1] can be obtained by results that are already 
available in the literature, as indicated in the theorems presented in this paper. Most results in the 
literature concern about the convergence of [L/M\ Pade approximants to a Stieltjes function as L 
and/or M — > oo. Naturally, they exclude those functions whose spectral function \i takes finitely 
many different values due to the fact that these functions are rational functions and [L/M] Pade 
approximants for a rational function is the function itself when L and M are large enough; hence 
convergence is guaranteed. The motivation for proving the main result presented in Theorem 3. 1 is 
due to the concerns on its application on inverse homogenization as described in the introduction. 
It's well known that for any rank n laminate microstructure, F(s) has finite number of poles, i.e. 
its spectral measure takes finite number of different values. At the same time, it is also known 
that the checkerboard microstructure has a spectral function which takes infinitely many different 
values and has a jump at z = 0. In other words, given data on effective dielectric constant of a 
composite, we only know that the corresponding F(s) function can be represented as a Stieltjes 
integral but we do not know a priori how many different values the spectral function takes because 
the microstructure itself is the unknown for the inverse problem. From the inverse problem point 
of view, the existence of the nonstandard Pade approximants is as important as the convergence 
property. As mentioned in the introduction, in Q, the spectral measure was reconstructed with 
very good precision by using data of effective dielectric constants at different frequencies. It is 
done by recovering the coefficients bo, 62, . . . , &m> a<o, ■ ■ ■ , o-M by reformulating the problem as 
minimization problem with constraints on the location of poles s n and residues A n in Q. The 
results in this paper also justify these constraints except Y^ n A n < 1 by identifying A n with — ^ 



and s n with —l/p n in (28 ). The constraint Yl n < 1 is justified by noticing that ^ is the 



Pi 



constant term of the Maclaurin series of [M/M — 1] j in (26 1 and it has to be equal to the constant 
term of the Maclaurin series of /, which is /j® in ([6]). fXo < 1 because it is the volume fraction of 
one of the constituent materials. 
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